Coherent Quasiclassical Dynamics of a Persistent Current Qubit 
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A new regime of coherent quantum dynamics of a qubit is realized at low driving frequencies in the 
strong driving limit. Coherent transitions between qubit states occur via the Landau- Zener process 
when the system is swept through an energy-level avoided crossing. The quantum interference 
mediated by repeated transitions gives rise to an oscillatory dependence of the qubit population on 
the driving field amplitude and flux detuning. These interference fringes, which at high frequencies 
consist of individual multiphoton resonances, persist even for driving frequencies smaller than the 
decoherence rate, where individual resonances are no longer distinguishable. A theoretical model 
that incorporates dephasing agrees well with the observations. 
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Macroscopic quantum systems coherently driven by 
external fields provide new insights into the fundamen- 
tals of quantum mechanics and hold_promise for applica- 
tions such as quantum computing Ij. Superconducting 
Josephson devices are model quantum systems that can 
be manipulated by RF driving fields , and recent years 
have seen rapid progress in the understanding of their 
quantum dynamics ^, ^ H, H 0, H .10]. Quantum 
coherence of these systems can be probed by temporal 
Rabi oscillations 0, E S S, H 0- There, the driving- 
field frequency u equals the energy level separation AE^ 
and the population of the two levels oscillates at a fre- 
quency cjR much smaller than AE. In the weak driving 
limit, hcjB. ~ A <C AE = hu^ where A is the driving 
amplitude parameterized in units of energy. 

Coherent quantum dynamics can also be investigated 
at driving frequencies much less than AE^ and at strong 
driving amplitude A ^ AE hv. In this case, the 
transitions occur via the Landau-Zener (LZ) process at 
a level crossing 0, ^3- Acting as a coherent beam- 
splitter, LZ transitions create a quantum superposition 
of the ground and excited states and, upon repetition, 
induce quantum mechanical interference. The latter 
leads to Stueckelberg-type oscillations 0, Q in anal- 
ogy to a Mach-Zehnder (MZ) interferometer 0, H^ . 
These oscillations are also related to photoassisted trans- 
port 0, 18, 19] and Rabi oscillations observed in the 
multiphoton regime 0, ji^. MZ-type interference is a 
unique signature of temporal coherence complementary 
to Rabi oscillations, with the time between sequential LZ 
transitions clocking the dynamics similarly to Rabi pulse 
width. 

In this Letter, we report a new quasiclassical regime 
which exhibits coherence even at driving frequencies low 
compared to dephasing rate, vT2 ^ 1 21]. This occurs 
because the interval between consecutive LZ transitions, 
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FIG. 1: (a) Schematic of the PC qubit surrounded by a 
DC SQUID readout. DC and RF fields control the state of 
the qubit. (b) An RF pulse of duration At ^ drives 
the qubit, and its state is inferred from the voltage Vsq 
across the SQUID pulsed with current /sq. (c) The qubit 
experiences two Landau-Zener transitions over a single RF 
period, accumulating a relative phase AO12 between them, 
(d) The resulting interference fringes in qubit population for 
u = 270 and 90 MHz, and V^ms = 240 and 171 mV respec- 
tively (vertical lines on FigEJ. 

relevant for MZ interference, is only a fraction of the driv- 
ing field period. We investigate the crossover between the 
multiphoton and quasiclassical regimes, demonstrating 
that coherent MZ-type interference fringes in the qubit 
population persist for frequencies UT2 ^ 1 even though 
individual multiphoton resonances can no longer be re- 
solved. This behavior should be contrasted with Rabi 
oscillations, where at low driving frequency, UT2 ^ 1, 
there is no signature of coherence. The crossover be- 
tween the two regimes, h'T2 ~ 1, is also influenced by 
inhomogeneous broadening, as discussed below. 

In our experiment we utilize a persistent-current (PC) 
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qubit l23|: a superconducting loop interrupted by three 
Josephson junctions (JJ), one of which has a reduced 
cross-sectional area (Fig^). A time-dependent mag- 
netic flux f{t) = /^^ + /^^ controls the qubit. For 
f{t) ~ ^0/2, the qubit exhibits a double- well poten- 
tial proflle with individual wells representing diabatic 
circulating-current states, |0) and |1), with energies ±e oc 
±(5/, where 6f = f^^ — <l>o/2 is the flux detuning. 
These states are coupled with a tunneling energy A. 
The driving and readout pulse sequence is illustrated in 
FigH^b). Qubit transitions are driven by a microwave 
flux /^^ (X Acos27rut^ with A, parameterized in units 
of energy, proportional to the microwave source voltage 
Vrms' The qubit state is read out with a DC SQUID, 
whose switching current Isw depends on the flux gen- 
erated by the qubit and, thereby, the qubit circulating- 
current state. The device was fabricated at Lincoln Lab- 
oratory using a fully-planarized niobium trilayer process 
and optical lithography. The device has a critical current 
density Jc ^ 160A/cm^, and the characteristic Joseph- 
son and charging energies are Ej ^ (27r/i)300 GHz and 
Eq ~ (27rfi)0.65 GHz respectively. The ratio of the qubit 
JJ areas is a ^ 0.84, and A ^ (27rfi)10 MHz. The ex- 
periments were performed in a dilution refrigerator at a 
base temperature of 20 mK. The device was magnetically 
shielded, and all electrical lines were carefully flltered and 
attenuated to reduce noise (see Ref. details). 

The qubit dynamics in the strongly driven limit is in- 
fluenced by quantum interference at sequential LZ transi- 
tions. As illustrated in FigQ^c), the qubit is initially pre- 
pared in the ground state at flux detuning Sf, and, after a 
flrst LZ transition at time ti, it is in a coherent superposi- 
tion of the two diabatic states. For times ti < t < ^2, the 
superposition state accumulates a relative phase A^i2, 
which mediates the quantum interference at the second 
LZ transition at time ^2- The sequence of two LZ transi- 
tions, repeated many times during the RF pulse, is anal- 
ogous to a cascade of MZ interferometers. One expects 
MZ-type interference fringes in the qubit population due 
to changes in A^i2 associated with changes in Vrms and 
6f, which are indeed observed (Fignjl). 

Figure [2l presents the measured qubit population of 
state |1) (color scale) as a function of Vrms and 6f for 
high- and low-frequency driving, u = 270 and 90 MHz 
respectively. Population transfer due to qubit driving 
appears at Vrms exceeding a threshold value which varies 
linearly with |^/| and symmetrically about the qubit 
step. For high-frequency driving, h'T2 ^ 1, the individ- 
ual multiphoton resonances are distinguishable and form 
a "Bessel ladder" 15j(Fig. [2^). The population of state 
|1) for the n*^-photon resonance follows a Bessel- function 
dependence, J^{A/hiy). The range of Sf in Fig. 2a ac- 
commodates photon transitions with n = 1 — 45, which 
together deflne coherent MZ interference-fringe bands of 
discrete resonances. 

In contrast, for low- frequency driving, UT2 ^ 1, the in- 
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FIG. 2: Measured qubit population at strong driving in two 
regimes, (a) v = 270 MHz. Multiphoton resonances of order 
up to n = 45 can be discerned (z^T2 > 1). (b) u = 90 MHz. 
Individual resonances are no longer distinguishable {iyT2 ^ 
1), but coherent interference is still observed. Vertical lines 
indicate the scans displayed in Fig^. A pulse of duration 
At = 3 /iS was used in both cases. 

dividual photon resonances are no longer distinguishable 
because the resonance widths exceed the resonance spac- 
ing (Fig. I^b). Nonetheless, the MZ interference- fringe 
bands, a signature of coherence in the strongly driven 
regime, indicate that the coherent interference mediat- 
ing the population transfer persists. 

To understand these results we consider a driven qubit 
subject to the effects of decoherence {h = 1): 

^ 2 V A -h{t) J ■ 

Here, h{t) = e + 5e{t) + Acos27:vt is the energy detun- 
ing from an avoided crossing modulated by the driving 
fleld in the presence of classical noise 5e(t). By a gauge 
transformation, the Hamiltonian is brought to the form 

where 4){t) = h{T)dT. Perturbation theory gives the 
rate of LZ transitions between the states |0) and |1): 

/•*' 

W= hm \At,t'?/5t, At,t'= / i^{t')dt', (3) 

where St = t' — t > 0^ and the limit physically means 
that 5t is large compared to T2. 

For the perturbation approach to be valid, the change 
of qubit population must be slow on the scale of T2 . This 
condition can be written down as W <^T2 = I/T2. We 
stress that this inequality does not imply that the effect 
of driving the qubit is weak. The rate W can still be 
large compared to the inelastic relaxation rate Fi, leading 
to the strong deviation of population from equilibrium 
observed in our experiment. 
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To evaluate VF, we write the expression in Eq.(|3|) as 
W=^J (A(t + r)A*(t))^^r, A{t) = Ae-'^^^'l (4) 
By introducing Bessel functions in the Fourier series of 

g-i(A/a;)sinu;t^ ^ ^ 27ri/, WC haVC 



-i(f)(t) —iet — i5c/){t) 



_A_A 
CO hv 



We average over 5(j)(t) with the help of the white noise 
model (e^^^(^)-^^^(^')) = e-r2|t-t'l^ and integrate in ^ 
as / e-^(^-^^)^-r2|r|^^ = 2r2/((e - unf + Vl) to obtain 



W{e,A) 



(e - wn)2 + Tl ■ 



(5) 



For large n, Bessel functions can be expressed through 
the Airy function A\{u) = ^ cos{uy + \y^)dy as 
Jn{x) = aAi{a{n — x))^ a = Using the iden- 

tity cotz = — 7rn)~^ we approximate Eq. (0) as 
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lmcot(-{e-iT2))Ai^(-{e-A)). (6) 



There are two main regimes exhibited by this expression: 
(i) ^ and (ii) u ^ In case (i), we have a sum 
of non-overlapping resonances. For each value of e, the 
sum is dominated by the term with n the nearest integer 
to e/cj, giving rise to resonances of strength J^(^), the 
Bessel ladder of Ref. . 

In contrast, in case (ii), the peaks in Eq. ((SJ are over- 
lapping. Setting cot = z in Eq.© 21], we obtain 



W{e,A) 
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Ai^{a{e- A)/uo). 



(7) 



The effect of r2 on the Airy function oscillation is small 
at r2 ^ (2i\la)v. Since a ~ 0.3 for ejlnv ^ 50, this 
condition is compatible with v ^^2- Eq.© can also be 
obtained by considering just two subsequent passages of 
a level crossing at a short time separation 1^2 ^ 
and ignoring the periodicity of the driving. 

Since Ai(ii < 0) oscillates as 7r"^/^|^i|"^/^ cos(||^ip/^ - 
^), while k\{u > 0) decays exponentially, Eq.© implies 
that the transitions occur only for A ^ e, with a rate 
which oscillates as a function of A — e. The oscillations 
are the same for both integer and noninteger e/hu^ con- 
firming that, while the resonances merge into a contin- 
uous band, the interference fringes persist at ^ r2, in 
agreement with our observations. 

To describe the population dynamics in the presence 
of driving, we employ a rate equation approach, in which 
the qubit level occupations obey Pi = QijPj , where 
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-911 = W + El, gio = -goo = W^ T[. (8) 



Here, Fi = 1/Ti, F; = Fie"^^ are the down and up 
relaxation rates. The magnetization of the stationary 
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FIG. 3: (a) Time evolution of excited state population = 
90 MHz, Vrms = 171 mV) obtained by varying the pulse width 
At. (b) The characteristic rate F as a function of flux detuning 
6f obtained by fitting to the exponential time dependence 
[Eq.Q] (inset shows examples of fits for the points I, H, HI 
and IV). 

state is rus = po - Pi = (Ei - T[)/{2W + Fi + F^), 
which gives the equilibrium value mo = tanh ^/3e at weak 
excitation, and <C mo at high excitation. 

To validate this model, we investigate the interference 
fringes in the excited state population as a function of 
RE pulse length At (EiglHt)- The rate equation predicts 
an exponential time dependence for the magnetization. 



m{At) — ms-\-{mo — ms)e' 



-FAt 



F = 2W+Fi+F;. (9) 



By fitting exponentials [Eq.®] to the qubit population 
at each flux detuning, we find the rate F which character- 
izes how fast the stationary state is approached (FiglHt)). 
Since our Ti ^ 20 /is is much longer than the ob- 
served transition time, we have F ^ 2W . Comparing the 
extracted F with Eq.Q), we obtain A/27r/i = 13 MHz. 

In Fig. ^ we compare the experimental characteris- 
tic rate and qubit population with those predicted by 
the model [Eqs.Q,®]. From the best fits we obtain 
F2/27r = 12 - 18 MHz (T2 ^ 9 - 13ns), consistent with 
the transition between the multiphoton and quasiclassi- 
cal regimes of Fig. El 270 MHz > F2 = I/T2 > 90 MHz. 

Inhomogeneous broadening is incorporated into the 
model by assuming a Gaussian broadening mechanism 
with standard deviation cr/27r = 40 - 45 MHz. The 
resulting power-broadened linewidth is approximately 
150 MHz, consistent with the linewidth observed in 
Fig. I2t- Best fits in Fig. ^ are obtained with slightly 
different values of F2 and a within the ranges above. By 
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FIG. 4: Comparison of experiment (blue) and theory (red), 
(a) The transition rate from the right half of Fig|2l3 fitted 
with r defined by Eqs.(|5|,®- (b) State |1) occupation taken 
from FigEk, compared to the model, Eq.Q. 

using the fit parameters for the 3/is magnetization curve, 
we can calculate the qubit population in the multipho- 
ton (Fig. [5^) and quasiclassical (Fig. [SJd) regimes as a 
function of 6f and Kms- 
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FIG. 5: Simulation of qubit population using model parame- 
ters extracted from data, (a) u = 270 MHz. (b) u = 90 MHz. 

In conclusion, we have observed quantum coherent 
qubit dynamics at strong driving for frequencies smaller 
than the dephasing rate. In this limit, well-resolved mul- 
tiphoton transitions merge into a continuous band, while 
the Mach-Zehnder-like coherent interference pattern per- 
sists. A simple model of a driven two- level system sub- 
ject to decoherence is in remarkable agreement with the 



observed interference patterns. 
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